Observables for bound orbital motion in axially symmetric space-times 
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The periastron shift and the Lense-Thirring effect of bound orbital motion in a general axially 
symmetric space-time given by Plebahski and Demiahski are analyzed. We also define a measure for 
the conicity of the orbit and give analytic expressions for the observables in terms of hyperelliptic 
integrals and Lauricella's Fd function. For an interpretation of these analytical expressions, we 
perform a post-Schwarzschild and a post-Newton expansion of these quantities. This clearly shows 
the influence of the different space-time parameters on the considered observables and allows to 
characterize Kerr, Taub-NUT, Schwarzschild-de Sitter, or other space-times. 



I. INTRODUCTION 

Axially symmetric vacuum solutions of the Einstein 
field equations are used to describe a wide range of black 
holes appearing in the universe. The most prominent 
of these solutions is the Kerr space-time [1], which de- 
scribes an axially symmetric rotating black hole. Gener- 
alizations of these are the charged Kerr-Newman black 
holes, the Kerr-de Sitter space-time which incorporates 
a non- vanishing cosmological constant, the Kerr- Taub- 
NUT space-time which includes the NUT-charge, or some 
combination of these. All these space-times can be seen 
as special cases of the family of electrovac space-times of 
Petrov type D given by the Plebanski-Demiahski class of 
solutions [2]. They are characterized by seven parame- 
ters: the mass, the rotation around the symmetry axis, 
the electric and magnetic charge, the NUT charge, the 
cosmological constant, and the acceleration of the gravi- 
tating source. 

Shortly after the introduction of general relativity, Ein- 
stein found that for every revolution of a planetary orbit 
the point of the shortest distance between particle and 
central object, the periapsis, is shifted in the direction 
of rotation of the test particle [3(. Indeed, the expla- 
nation of the anomalous shift of Mercury's perihelion, 
together with the observation of light deflection, consti- 
tuted the breakthrough of general relativity. Already in 
1918, Lense and Thirring found that bound orbital mo- 
tion around an axially symmetric rotating black hole is 
perturbed not only compared to the Newtonian case but 
also to the non-rotating case [4J . The rotation of the cen- 
tral object causes an additional shift of the periapsis and 
a precession of the orbital plane. For the more exotic 
Taub-NUT space-time Misner and Taub Q showed that 
the geodesic motion takes place on a cone which, if slit 
open and flattened, defines the orbital plane as it would 
be for vanishing NUT charge. These three examples show 
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that the parameters of the space-time affect the observ- 
ables of bound orbital motion and that in turn these ob- 
servables may be used to characterize a space-time. In 
this paper we will investigate these orbital effects starting 
from the general Plebanski-Demiahski space-time with 
vanishing acceleration of the gravitating source. 

The only way to get access to these characteristics of 
black holes is through the orbits of particles and light 
around these black holes. They are given as solutions of 
the geodesic equation describing the motion of test parti- 
cles and light rays in a given space-time. For a thorough 
investigation of the physical properties of orbits analyt- 
ical solutions of the geodesic equation are most useful. 
Hagihara [6[ was the first to find an analytical solution 
of the geodesic equation in Schwarzschild space-times us- 
ing elliptic functions. He also classified all possible types 
of orbits in terms of the energy and the angular momen- 
tum of the test particle. The complete set of solutions 
for charged particle motion in Reissner- Nordstrom space- 
times have been presented only recently by Grunau and 
Kagramanova [7| . Shortly after the discovery of the Kerr 
solution in 1963 [l|, a number of authors studied the 
geodesic motion in this space-time. Their results were 
reviewed and extended by Chandrasekhar 8]. A treatise 
of the geodesies in Kerr space-time of the same com- 
pleteness as 6] was given only recently by Slezakova Q 
and for the motion of charged particles in Kerr-Newman 
space-times by Xu [10]. Geodesies in even more com- 
plicated space-times like the one incorporated in the 
Plebanski-Demiahski space-time are beyond the methods 
introduced by Hagihara. However, in 2008 a method has 
been found to analytically integrate the geodesic equa- 
tion in Schwarzschild-de Sitter space-times using the the- 
ory of hyperelliptic functions [FL], [l2[ and to classify all 
types of possible orbits in terms of energy and angular 
momentum of the particles as well as the cosmological 
constant. This approach was then used to find solutions 
of the geodesic equations in general Plebanski-Demiahski 
space-times with vanishing acceleration of the gravitating 
source [l3| and to some special cases [14], [l5| . 

Bound stable orbits far away from the black hole are 



not very different from Newtonian orbits but in the vicin- 
ity of the black holes they exhibit some peculiar features. 
These are caused by mismatches between the periodici- 
ties of the radial, polar, and azimuthal motion, which 
are all equal to 2-7T in the Keplerian case but may differ 
greatly near black holes. A systematic approach to han- 
dle these different periodicities for Kerr space-time was 
worked out by Schmidt [16] and continued by Drasco and 
Hughes [17| using a technique for decoupling the radial 
and polar motion suggested by Mino [18]. They showed 
that arbitrary functions of Kerr black hole orbits can be 
described in the frequency domain and how these ex- 
pansions may be computed explicitly. Later, Fujita and 
Hikida [19| derived analytical expressions for the frequen- 
cies of radial, polar, and azimuthal motion in Kerr space- 
time and also for bound timelike orbits in Kerr space- 
time. 

In this paper, we will find analytic expressions of the 
fundamental frequencies defined by Schmidt [16j , Drasco 
and Hughes [l7|, as well as Fujita and Hikida [19j for 
the general Plebahski-Demiahski space-time with vanish- 
ing acceleration of the gravitating source. We will then 
investigate how these frequencies are influenced by the 
parameters of the black hole. To this end, we expand 
the analytic expression in terms of the parameters in a 
Taylor series up to first order. This clearly shows the in- 
fluence of different parameters on the orbits of the black 
hole and leads to conclusions which kinds of orbits may 
belong to a given black hole. 

The paper is organized as follows: First, we review 
the equations of motion in Plebahski-Demiahski space- 
times decoupled by the method introduced by Mino. 
These equations are then used in the third section to 
define the fundamental frequencies and observables for 
Plebahski-Demiahski space-time following closely the ar- 
guments of Fujita and Hikida. All necessary quantities in 
Schwarzschild space-time as well as their post-Newtonian 
expansions are computed in section four. These quan- 
tities will serve as a reference for the comparison with 
more complicated space-times later on. Section five is 
the most technical one, where the linear correction to 
the fundamental frequencies are computed. These cor- 
rections are given in terms of elementary function or in 
terms of complete Jacobian elliptic functions. In section 
six we use this results to compute the post-Schwarzschild 
and post-Newtonian corrections to the periastron shift, 
the Lense-Thirring effect, and the conicity and compare 
them to earlier results. A discussion and outlook closes 
the paper. 



II. GEODESIC MOTION IN 
PLEBANSKI-DEMIANSKI SPACE-TIME 

The axially symmetric Plebahski-Demiahski space- 
times are characterized by the seven parameters mass, 
rotation, acceleration, cosmological constant, NUT- 
parameter, electric and magnetic charge. It can be 



shown that the Hamilton- Jacobi equation for these space- 
times is separable and the geodesic equation integrable, 
if and only if the acceleration vanishes or null geodesies 
are considered. In this paper we will consider the 
case of geodesic motion of massive test-particles in a 
Plebahski-Demiahski space-time with vanishing acceler- 
ation. We will also assume that the test-particles are 
neutral, i.e. without electric or magnetic charge. 

The six-parameter Plebahski-Demiahski space-times 
considered here are then given by the metric @, [20] [21] 
(we use units where c = 1 = G) 



ds 2 /M 2 
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where p 2 = r 2 + (n — a cos 6) , A — a sin 2 6 + 2ncos9, 
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Here a, the angular momentum per mass of the gravitat- 
ing source, the NUT parameter n, the electric charge Q c , 
and the magnetic charge Q m as well as the coordinates 
r and t are normalized with respect to M, where M is 
the mass of the gravitating object. The dimensionless 
parameter A denotes the cosmological constant divided 
by three and normalized by multiplication with M 2 . 

The equations of motion for massive test-particles in 
these space-times are given by 



|1) -/Vr A,(,- ff>: :/?[,-) (4) 



Q{9) 2 
sin 2 9 



A e (C - (n - acos6>) 2 ) - — '-+- =: 6(0) 
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dt r 2 + a 2 + n 2 , . a sin 2 9 + 2ncos9 „,„. 

IT = * p ( r ) + ; 5 0(6) 

d\ A r w A e sin 2 6> V ' 



where 



--:T(r,6) 

P(r) = (r 2 + a 2 + n 2 )E - aL , 
0(6) = L- (a sin 2 6 + 2n cos 6)E . 



(7) 

(8) 
(9) 



The constants of motion E, L, and C have the meaning 
of energy, angular momentum in direction of the symme- 
try axes, and Carter constant, each per unit mass. The 
two constants L and C are additionally normalized by 
division by M and M 2 , respectively, such that they are 
dimensionless. The afiine parameter A is the Mino time 



18], also normalized to M, dX = dr/(Mp 2 ) with the 
eigentime r. Note that R depends quadratically on n, 
Q c , and Q m and that 9 does not depend on Q e and Q m . 
The equations (H)-® can be solved analytically [1J, 
[22|, [23|. However, in this paper we are interested in the 
periods of the motion, which can be used to define ob- 
servables related to geodesic motion in these space-times. 



III. OBSERVABLES FOR BOUND ORBITS 

In Newtonian gravity bound geodesic motion is de- 
scribed by a fixed ellipse defining the orbital plane. This 
is no longer true in the framework of general relativity. 
In the spherically symmetric Schwarzschild space-time, 
the orbital plane remains fixed while the ellipse precesses 
resulting in a periastron shift. In the axially symmetric 
Kerr space-time additionally the orbital plane itself pre- 
cesses what is known as the Lense-Thirring effect [24]. In 
the more exotic Taub-NUT space-time the geodesic mo- 
tion of test particles does not lie in a plane at all but on 
a cone, see e.g. [251 ] . 

The precession of the orbital ellipse and the orbital 
plane is induced by mismatches of the periods of the mo- 
tion in the r and coordinates compared to the average 
secular increase of the angle tp about the symmetry axes. 
These effects were discussed in the framework of a Kerr 
space-time by Schmidt [lj| , Drasco and Hughes [17| , and 
Fujita and Hikida [19j . In the following sections, we will 
use the procedure in [19J ] to analyze the first order cor- 
rections to the periastron shift and the Lense-Thirring 
effect due to the parameters a, n, and A. (As Q e and 
Q m appear only quadratically, there are no linear effects 
due to electric or magnetic charge on neutral test par- 
ticles.) Below we review their line of argument for the 
convenience of the reader. In addition, we will character- 
ize the deviation from an orbital plane, the conicity, due 
to the parameter n. 

For bound orbits, the radial and polar components r 
and vary between a minimal and maximal value given 
by the turning points ^ — and J| = 0. The periods 
A r of r and Ag of 9 with respect to the Mino time A are 
then defined by a revolution from maximum to minimum 
and back to the maximal value. This means that A r and 
Ag are defined by the smallest non-zero real value with 
r(X + Ar) = r(A) and 9(X + Ag) = 0(A) giving 



Ar = 

A s = 



dr 



= 2 



dr 
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d9_ _ f ""** d9 



(10) 

(11) 



where r p is the periapsis and r a the apoapsis. In general 
these are hyperelliptic integrals, which are a generaliza- 
tion of elliptic integrals. (The closed integration path o r 
(ag) refers to the integration on the Riemann surface of 
the algebraic curve defined by y 2 — R(r) (y 2 = 0(cos0)). 
It runs around the branch cut connecting r p and r a 



(cos0 m i„ and cos9 max ). On the Riemann surface the 
two branches of the square root are glued to one analytic 
function.) From these two periods conjugate fundamen- 
tal frequencies can be defined by 



T. r := signL — 



2tt 
Tg := sigmL — . (12) 

Jyg 



The sign of L is included in this definition to indicate the 
direction in which the particle travels around the gravi- 
tating object relative to a > 0. This means that L > 
corresponds to a prograde and L < to a retrograde 
orbit. 

For vanishing n, Eq. ([5]) is symmetric with respect 
to the equatorial plane = j. However, in the more 
general space-time considered here this is not true, and 
the deviation from this symmetry can be measured by 
the difference between \{9 min + 9 max ) and ?. In the 
case of a Taub-NUT space-time, where all parameters 
except mass and NUT parameter vanish, this implies 
that the particle moves on an orbital cone rather than 
an orbital plane pj. Later we will see that this phe- 
nomenon also appears in the weak field limit. For these 
reasons, we will refer to this quantity as the conic- 
ity, A conicity := 7r - (0 m , ; „ + 9 max ). This means that 
A con i C ity > corresponds to a cone opened in northern 
direction and A con i C j tv < to a cone opened in southern 
direction. 

The nature of the equations (JU) and (J7J is somewhat 
different from (fJJ and ((5|) as they cannot be solved by 
periodic functions. They depend on both r and but 
can be separated in an r-dependent and a 0-dependent 
part 

$(r,0)=:$ r (r)+$ (0), (13) 

T(r,9)=:T r (r)+T e (9). (14) 

The solutions <p(\) and t(X) of equations ([5]) and ([7]) can 
be written as an averaged part linear in A plus perturba- 
tions in r and 

<p(\) = (fc(r, 0)) aA + *£. c (r) + $L(») (15) 

t(X) = (T(r, 9)) x \ + T r osc {r) + T e osc {9) (16) 

where 

( ' >A:= a ^ 2fA l X) f^^ (1?) 

(A 2 -Ai)-s-oo Z(A 2 — M) JXi 

is an infinite time average with respect to A, and X^ sc (r) 
and Aq SC (0) (with X = $ or X — T) represent oscillatory 
deviations from this average. They are defined by 



Ksc{r)= X r (r)dX-(X r (r)) x X. 



X e osc (9)= Xg(9)dX-(Xg(9)) X X 



(18) 
(19) 



and have periods A r and Ag. Therefore, the average sec- 
ular increase of tp and t with respect to A is given by 

r:={T(r,e))x = {T r (r))x + {T d (9))x, (20) 

T v := ($(r, 0))a - (* r (r))A + (**(0)>a , (21) 



As X r and Xg (again with X = $ or X = T) are pe- 
riodic functions with respect to A because r and 9 are 
periodic functions, their integrals § x 2 X r (r(X))d\ and 

f, 2 Xg(9(X))dX in the definition of the infinite time av- 
erage can be reduced to an integral over one period A r 
or Ag, respectively. This yields 
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' ArJ rp y/Wjr) 
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Ae 



<S>g(0)dO 
Tg(9)d0 
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(22) 
(23) 



The corresponding frequencies with respect to the coor- 
dinate time t are then given by 






r 



n ( , 



r 



a 



r 



(24) 



These are frequencies as seen by an observer at infinity. 
In the limit of weak gravitational fields the mismatch of 
the frequency of the ip- and r-motion, fip = O v — Q r can 
be interpreted as the precession of the orbital ellipse, and 
the mismatch of the frequency of the if- and 0-motion, 
£1lt = Qip — £^g as a precession of the orbital plane. In 
strong gravitational fields the orbits are in general of such 
irregularity that orbital planes or ellipses can no longer 
be identified. 



(E 2 - l)r(r - r 02 )(r - r 03 ){r - r 04 ) yields 



E z - 1 = 



C = 



ro2 



ro2 + r 03 + r 04 
f027-03n>4 

7"02 + ^03 + r 4 ' 
-2r 03 r 4 

2r 03 + 2r 04 - r 03 r 4 



(27) 
(28) 
(29) 



A non-polar orbit requires C > L 2 and lies in an orbital 
plane with inclination arcsin ±y= (or -| — arcsin ±-X if mea- 
sured from the equatorial plane). Therefore, the theta 
motion of the reference orbit is symmetric with respect 
to the equatorial plane and confined to [0oij$O2] with 
0oi = arcsin -LJ= g [0, §] and 6> 02 = 7r-arcsin ^ e [f i ""]■ 
In particular, the orbit lies in the equatorial plane if 
C = L 2 . In terms of the inclination and the turning 
points the constant of motion L is given by 



L = ±\/Csin6> i =± 



^02 + r 3 + r 04 



sin W01 



(30) 



In the following we will calculate the expressions T r , Tg, 
T^, and T for the reference orbit. 



A. Frequency of r 



IV. REFERENCE ORBIT 

In the next section we will analyze the influence of the 
parameters a, A, and n on the observables of a bound 
reference orbit, which is neither circular nor polar. (Po- 
lar orbits may be considered as a special case but we 
leave this out here. For a discussion of these orbits in 
Kerr space-time see |2g.) This will be done by post- 
Schwarzschild and post-Newton expansions, where we as- 
sume the constants of motions to be fixed. But first we 
will introduce and characterize the reference orbit. 

If all space-time parameters except the mass vanish, 
the Plebahski-Demiaiiski space-time is identically to the 
Schwarzschild space-time and the functions R(r) and 
6(0) of Eqs. and (© reduce to 



R (r) := (E 2 - l)r 4 + 2r 3 - Cr 2 + 2Cr , 
e (6):=C- L " 



sin 



(25) 
(26) 



A necessary condition for the existence of a non-circular 
bound orbit in a Schwarzschild space-time is that Ro(r) 
has four real zeros, — rr,i < rr,2 < ro3 < rr,4 < 00, 
where Ro(r) > for r 03 < r < r 4, cf. [23j. Therefore, 
our bound reference orbit has turning points ro3 and tqa- 
Note that instead of using E 2 and C all formulas can also 
be expressed in terms of the turning points: A compar- 
ison of coefficients in (E 2 - l)r 4 + 2r 3 - Cr 2 + 2Cr = 



Let us first calculate the r period for our reference 
orbit. For an orbit bound between ro3 and ro4 we get 



lj-,0 



dr 



VR^r) 



(31) 



This is a complete elliptic integral of first kind which can 
be easily transformed to the Legendre form giving 



^r\0 



4K(k) 



i/(l - E 2 )r Q3 (r 4 - r 02 ) ' 



(32) 



with k 2 = r 2(r i ~ r 03 ) / (r 03(^4 ~ r 02 )). For general 
information on the complete elliptic integral of first kind 



Kih) = Jo7n 



dt 



^{l-t 2 )(l-kH 2 ) 



(33) 



see e.g. [27J, for fast numerical computation see e.g. |28| . 
In computer algebra systems like Mathematica or Maple 
the complete elliptic integrals are usually implemented 
and, therefore, A r $ can be computed easily. The conju- 
gate fundamental frequency T r! o is then given by 



T r .o — 



2tt 



AK{k) 



sign! i/(l - E 2 )r 3(r 04: 



ro2) 



(34) 



B. Frequency of i 



D. Frequency of t 



The 9 period of the reference orbit is given by 

f e ° 2 dd 
Ag. := 2 



e 01 VQoJo) 

?03 d£ 



?02 ^c{i-e)-L^ 



(35) 



where we substituted £ = cos and cos #oi = \ 1 



£ 3, COS 6*02 =6 



L-> 



02 



-£o3- This can be solved by 



. 27T 27r % /r 02 + r 03 + r Qi 
A e,o - —/= — . 



from which we infer 



■0,0 



signZ 



2tt 



Note that for the case 0qi 



= sisiiLVC. 



% 2 it is 0(A) 



(36) 



(37) 



and, 



therefore, Ag$ is undefined. However, we can treat this 
as the limiting case C — > L 2 , which gives the same results 
for Agp and Tg t o as above. 



C. Frequency of ip 

For the reference orbit Eq. (J6j> simplifies to 



<kp L 

dX ~ sin 2 



(38) 



and T v can be calculated to 



T 



tp,0 



L 



2n signL 
Ab,q 



signL v C 



(39) 



where we used the substitution £ = cos as in sec. IIVBI 
Although Agp is not defined for 0(A) = § we get in this 
case y(A) = LA and, therefore, the same result T v ,o = 
L = signL VC. 



The expression for To is the most complicated in this 
section as it involves an elliptic integral of third kind. For 
the reference orbit Eq. simplifies to 



dt r 3 E 



which leads to 



r = 



d\ r-2 



l Edr 



(40) 



A rfi J ro2 (r-2)/RoW 
E 



(41) 



2K(k) 

2r 0i (3-2E 2 ) 
1-E 2 



- ro3r 04 K(k) + r 03 (r 4 - r a2 )E(k) 



U( ni ,k) + ^±-U(n 2 ,k)) 
r 4 - 2 



where E(k) and H(n,k) are the complete elliptic integrals 
of second and third kind, 



E(k) = I 1 V=£± 2)dt 

Jo v/(l-i 2 )(l 



U(n,k) 



V '(1 - i 2 )(l - k 2 t 2 ) ' 
dt 



(42) 
(43) 



/o (l-nt 2 )^(l-t 2 ){l~k 2 t 2 ) 

and the parameters m, n 2 are given by 

7-03 - fQ4 2(r 04 - r 03 ) 
m = < , ri2 = -, ^ < . (44) 



'03 



nD3 (^04- 2) 



E. Observables 

Let us now collect the results so far in this section. 
As expected for a — 0, the frequencies Tgo and T^o 
coincide and the Lense-Thirring effect vanishes 



^LT,0 



TfVo — Too 



0. 



(45) 



Likewise, for n = the particle moves on an orbital plane, 

Aco„icity,0 = 7T - (001 + <?0 2 ) = . (46) 

The periastron shift ilpg = (T^.o — ^r,o)Lo~ is given in 
terms of elliptic integrals 



^p.o 



signL 



2VCK(k) - ir^/r 03 (r 04 -r 02 )(l-£ 2 ) 



E 



(-r 03 r 0i K(k) + r Q3 (r Q4 - r 02 )E(k) + 2r "^ g f 2) n(n 1; k) + ^n(n 2 , fc) 



(47) 



The usual result for the periastron shift in terms of radi- the averaged X(ip) — y>T * we get the period of the 
ans or degrees can be found in the following way: using 



motion in terms of ip, r v :— r o A : ip \-t r(ipT 1 ), by 
observing that r v (tp + T v A r ) = r(pT^ + A r ) = r v ((p). 
From this period we get for the difference between the 
angle of the periapsis and 2n after one revolution 



Ap ; o = T Vi o^r,o — 27rsignL 



signL 



(48) 



4VCK(k) 



\/{^- E 2 )r m (r 4-ro 2 ) 



-2tt 



If this shift should be expressed in terms of radians per 
time, it can be referred to the time needed for a rev- 
olution from some fixed tp to <po again, i.e. the side- 
real period, or for a revolution from some fixed r to r 
again, i.e. the anomalistic period. The first choice cor- 
responds to the usual notion of the periastron shift in 
arcseconds per century whereas the second corresponds 
to the definition of $7p: the time elapsed for a revolu- 
tion from, say, periapsis to periapsis is given as the pe- 
riod of rt ■= r o A : t i-> r(tr _1 ), which is TA r . If 
Ap = T v A r — 27r signL is divided by this period we ob- 
tain n P = A P /(TA r ). 

Note that (|4"T1) and (|4*5| are exact and, therefore, more 
complicated than the post-Newtonian formula given 
e.g. in [29]. If we consider the weak field approxima- 
tion by assuming that the periapsis ro3 and the apoapsis 
ro4 become large, we recover the formula (51) of [29j with 
vanishing solar quadrupole momentum and PPN param- 
eters corresponding to general relativity, 



67r signL , r 



(49) 



(r a + r p ) is the semimajor axis and e 



where d = 

fj (?"a — fp) the eccentricity with the apoapsis r a = tqa and 
the periapsis r p = r 03 . Here we used (|2"T|) - (f2"9"|) to perform 
the series expansion. Note that the angular momentum 
of the test particle is usually chosen to be positive and, 
therefore, its sign does not appear in Eq. (51) of [29j. But 
we include it here as we choose the sign of L relative to 
the angular momentum of the gravitating source, which 
will be nonzero later on. The post-Newtonian expression 
for J7p is given by 



n P0 * 3signL Mi 

' dt(l-e 2 ) 



(50) 



which includes the perturbation of the time needed for 
one revolution from periapsis to periapsis. 

Let us test our results for the orbital motion of Mer- 
cury. From |30J we take for peri- and aphelion the values 

r min = (307500.7 ± 3.0) x 10~ 6 AU , (51) 

r max = (466696.6 ± 2.6) x 10~ 6 AU , (52) 

and use M — 1476.625 28 m for the mass of the Sun to 
determine r 03 = ^ff^ and r 4 = Ls ff L - With Eqs. ((27)) - 
(l29l and the results of this section we obtain 



n 



P,o 



3.252308 x 10~ 19 ± 5.5 x 10" 



21 



(53) 



In order to express this dimensionless quantity in terms 
of arcseconds per century we have to multiply it by -j^ 
where c = 299792458 ms _1 is the speed of light. Then 
we get 



f7 po _ = 42.980 48 ± (0.73 x 10~ 3 ) arcsec/cy, (54) 



in consistency with observations, cf. 31]. However, it 
is not the usual result in the sense that it describes the 
perihelion shift per revolution from periapsis to periapsis 
rather than from to 2-k. 

The usual result for the perihelion shift in terms of 
arcseconds per century can by found using Ap, 



A P o = 5.018 648 5 x 10~ 7 ± 3.7 x 10~ 12 . 



(55) 



In 30] a revolution of Mercury is given as 87.969257 days. 
We get with t^ := t o A : ip >— > T T~ ip for a revolution 
of 2?r 



27rr n M 



Ym 



TfVo c 



— =87.969 25 ±(8.5 x 10" 4 ) days, (56) 



which agrees to the given accuracy with observation. 
This yields for the perihelion shift 



Ap, 



42.980 48 ± (0.73 x 10~ 3 ) arcsec/cy , (57) 



in accord with observations. We see here that this 
value coincides with (|54[) within the given accuracy and, 
therefore, the two different definitions can not be distin- 
guished. 



V. FIRST ORDER CORRECTIONS 

In the following we will calculate the linear post- 
Schwarzschild corrections for all quantities used to de- 
fine the observables fip and J7lt as well as the conic- 
ity A con i C it y due to the parameters a, n, and A. As 
the parameters Q B and Q m appear only quadratically in 
Eqs. (|4]) to ([7]), we will not study them here. However, 
this would be a totally analogous procedure. We assume 
the constants of motion to be fixed but let the zeros Ti 
of R and 9i of vary. By this procedure we will reduce 
the hyperelliptic integrals appearing in the definitions of 
A r , Ag, Ttp, and T for the general Plebanski-Demianski 
space-time to elliptic integrals and elementary expres- 
sions. 

The observables flp and £7lt are defined through the 
frequencies T r , T#, T v , and T. Whereas T r and Tg 
depend only on the variable indicated in the index, T v 
and r can be separated in an r- and a ^-dependent part, 



Try* | ry 



ipO 



andr 



Tg , where 



T vr :- 



$ r (r)dr 

r7 rp V^W 



I, 



ipr 



T - e : ~ X 



r r := 



— $ 9 (fl)rffl = V 

,„ VW) ''' A o 

' a T r (r)dr _ _ J^ 
_2_ f 9 ™* T e (0)d# _ J tf 



(58) 



(59) 



VW) 



An 



In terms of these expressions the periastron shift and the 
Lense-Thirring effect can be restated as 



Q P = 



T v - T r I vr A 8 + I v gA r - signL 2irA e , . 
f :: - r a I r . ^ l,,l)l 

T^-T 



r 



J tr yl<? + I te A r 

I V rAg + I v9 A r - signL 2-kA t 

I tr Ag + ItgAr 



(61) 



A. Standard form of hyperelliptic integrals 



All integrals appearing here have the form 
^ f(x)dx 



VW 



(62) 



where Pq is a polynomial of degree 6 in x = r or 
x = £ = cos#, /(a;) is a rational function, and xi, X2 
are zeros of Pq. This type of integral is called hyperel- 
liptic integral. The functions P& and / as well as the 
zeros Xi, xi depend on the parameters a, n, and A. For 
a Taylor expansion of such an integral it is of advantage 
to first reduce it to a standard form similar to the elliptic 
integrals which appeared in the calculations for the ref- 
erence orbit. However, to our knowledge such a standard 
form does not exist in the literature. A straightforward 
generalization of the Legendre standard form of elliptic 
integrals of first kind can be obtained by an additional 
term (1 — k%t 2 ) under the square root but this yields again 
an elliptic integral as a substitution s — t 2 shows. A bet- 
ter choice is to generalize the Riemann form of elliptic 
integrals, 



dt 



/o y/t{l - t)(l - Wt) 
to the hyperelliptic form 

(At + B)dt 



y/t(l - t)(l - k 2 t)(l - k 2 t)(\ - k 2 t) ' 



(63) 



(64) 



The transformation of the form (J62I) to the form ([64)1 
depends on the range of integration [xi,X2]. In the case 
of Pq(x) = R(r), we have 4 real zeros r\ < r-i < r 3 < r^ 
and another two zeros tq , r§ which may be complex and 
which tend to infinity for vanishing A. Here x\ — r 3 and 
X2 = ^4 result in a transformation r = (At + B)^ 1 + n 
which yields 



dr 



VW? 



(At + B)dt 



D Jo yjt(l - i)(l - k'{t)(l - k 2 t)(l - k 2 t) 



(65) 



-*=lF D (±,&,2,0i 



n nF D f ^,PiA,rn 



where B — l/(f4 — ri) 



B r 



A = 



ri -r 3 



D = A 



(ri - n)(r 3 -n) 
{rz -ri)(r 4 -r 2 ) 



(r 4 -ri) 5 



k\ = 



-'■ A r , 




(66) 


(ri -r )(r 5 - 


-n) = 


:D r , (67) 


- k 2 




(68) 


— k 2 

— ■ K-2r ) 




(69) 


— k 2 




(70) 



(U -r 3 )(r 2 -Ti) _ 
(rz - n)(n -r 2 ) 

fc2 = ( r i-r 3 )(r 5 -n) 
(r-3 ~ n)(r 5 - r 4 ) 

k 2 = (r4-r 3 )(n - r Q ) 
(rz - ri)(n - r ) 



and Pi = (|, |, |), m = (fc^, fc| r , fe| r ) =: m 0r . For the 
case of Pg(x) = 6^(C) (see Eq. (|90p) there are two real 
zeros £2 < £3 and four other maybe complex zeros £0 , £1 , 
£4, £5, where £1, £4 tend to infinity for vanishing a and 
£0, £5 for vanishing a or A. With xi = £2 and X2 = £3 we 
get the same form ([6"5"j) but with _B = 1/(^3 — £1) =: B5, 



A = 



(6 - 6) 



«2-6)(&-Cl) 



^ = a 4 A (6_ei)(€4 



e 3 



(6 -a) 2 

, 2 _ fo-6)(c 4 -a) 

(6-ei)(e4-&) 

t2 _ (6-6)fe-a) 

/Co — 



A f 



■(£3-£o)fo-£5) 



.Dp 



A- 2 

ft l£ 



(6 

(£3 



■6)(6 
6)(£i 



(6 -6) (6 



and m = (fc 2 j, 



fr 2 t 2 ^ 



£3) 
Jo) 

£0) ■ 

=: m 0? 



A- 2 



A- 2 

K 3? i 



(71) 

(72) 
(73) 

(74) 

(75) 



which was also used in [26( for the calculation of the pe- 
riastron shift of equatorial orbits in Kerr-de Sitter space- 
times. As pointed out in [20|, the form (jMJ) can also 
be expressed in terms of Lauricella's hypergeometric Fn 
function, see also appendix [Al 



B. Linear correction to period of r 

In this section we will calculate the linear corrections 
to the reference orbit due to the parameters o, n, and A. 



The exact formula for the r period is given by 



we compare the coefficients of the equation 



Ar = 2 



dr 



3 VW) 



(A r t + B r )dt 



A- Jo jt{l-t)(\-kl r t){l-kl r m-klrt) 



2A r 7T /3 -• 

F D l-,^2,rn Qr 



2B r (\ r. 

~77T n \ 2 ,m ° r 



(76) 



All quantities indexed with an r are given as some com- 
bination of the zeros of R(r) and, therefore, depend on 
the parameters a, n, and A. This means that for a Tay- 
lor expansion of A r we need to explicitly know how the 
zeros Ti of R(r) depend on the parameters. To see this, 

I 



R(r) = A [r 



VA 



(r -n)(r -r 2 ) 



x (r - r 3 )(r- r 4 ) ( r 



(77) 



where fo = rovA, r§ = r^y/A do not have a singularity 
at A = 0. As the coefficient of r 5 vanishes in R(r) we 
immediately see that 



-(ra + r 3 + r 4 + n)A - (fo + h)VX = 



(78) 



and, thus, fo and fs have to be expanded in terms of vA. 
Therefore, we introduce I = yA as expansion parameter 
instead of A. Expanding the right hand side of (|TT|) in 
terms of a, n, and I gives a system of equations which 
can be solved for the derivatives of all r^ with respect to 
the parameters. 

The Taylor expansion of A r in p — (a, n, I) near p = 
then reads 



A r ^A r (p = 6) + Y, A-Ap = 0) 



! = 1 



(79) 






2£>J 



A r -F D f -,/3i,2, m 0r j + B r irF D f -,/3i,l,TOo, 



+ 



/d; 



Ar^rfc ( -,/3i, 2, m r] +B r ,iirF D f -, ft, 1, m , 






■ ; C1.2 



(kjr), 



^j^i 2k o r 



A- (f d (-J{,2,rh 0l 



-Fd ( -,/3i,2,m 0r j j +S r 7r f Fd f -,^,l,m 0r j - Fd f -,/3i,l,m , 



Pi, 



(80) 



p=0 



where we used ^|-Fd (a, /?, 7, af) = ^-(F D (a, P l ,j,x) — F D (a, 0,j, x)) which can be found in [33]. The symbol /3 J 

means that the jth component of /3 is increased by one (e.g. (3 2 = (|, I, 5)), and by Xj we denote the derivative of 
X with respect to p,. Note that 

<g P , forp = 0, (81) 



,2 __ ? '03 - fQ4 _ 2 



v 2r 



''03 



and, therefore, the Lauricella function Fd reduces to an elliptic integral in this case. Here r$i are again the zeros of 
i?o, i-e. the turning points of the reference orbit. 

For the linear correction with respect to a we obtain 



A r , a (0) = -4EL 



Ci(l - E 2 ) 2 (r 03 - r Q2 )K(k) + 2r 03 (4 - 3C + 3CE 2 )E(k) 
(1 - F 2 ) 3 v / r 03 (ro4 - r 02 ){l - E 2 )(r Q4 - r 02 )(r 03 - r 02 ) 2 (r 4 - r 03 ) 2 



where k 2 — ro2 )'' 04 _'' 03 ; as for the reference orbit and £1 is a non-symmetric function of rp 2 , tq 3 , and ro4, 

Ci = r Q2 r Q3 + r a2 r a4 - r 2 3 - r^ . 
The linear corrections with respect to n vanishes, 

A r ,„(0*)=0 



(82) 

(83) 
(84) 



and also the linear correction due to I = y/A, as expected. To determine the linear correction due to A we use 



— (a,n,A)| A=0 



Jl {a ' nJ) 2i 



(85) 



;=n 



In this way the linear correction due to A can be calculated to 



At, A (JS) 



\A-o3(ro4 -r 02 )(l-£ 2 )(l - E 2 ) 



^03?"04(7l - 2C72) 



(^02 - r 4)(r 03 - r Q4 ) 2 (r 02 - r 03 ) 



K(k) 



4r 03 C(-16C + 76CE 2 - 57C£ 4 - C 2 + C 2 E i - 48) 
(702 ~ r Q 4)(r Q2 - r 03 ) 2 (r 03 - r 04 ) 2 (l - E 2 ) 5 



E(k) 



6r 04 
(1-E 2 y 



-U( ni ,k) 



(86) 



where 71 and 72 are non-symmetric functions of ro2 , ?*03 , 
and ro4, 



71 = 2r 03 r 04 - 3r 03 r 04 - 3r 04 r 03 + 3r 03 r 02 

+ 3r^ 4 r 02 - 7-03 r 02 + 2 ''03?'04?'o2 ^ r 04 r 03f02 
~ r 04 r 02 - ^03 r 04?-02 , 

72 = r Q4: r Q 2 + r 03 r 2 - 2r 04 r 3 • 



(87) 



C. Linear correction to period of 6 

The linear correction to the 9 period can be determined 
analogously to the foregoing subsection. The exact for- 
mula for the period of the 9 motion is 



A f , 



1Ap 



+ 



/L\2 
2B f 



f d ^,Al,2,to 05 



/Dp 



kEd -,/3i,l,m 0€ 



Again, we have to determine the dependence of the zeros 
£i, i = 0, . . . , 5 on the parameters a, n and A. This time 
we use the ansatz 



I 

and £1 = a£i, £4 = a£ 4 , £ = ay/A£o, £5 = aVAt; 5 behave 
regular in the limit a = 0, A = 0. By comparing the 
coefficient of £ 5 it can be seen that £0 and £5 expand in 
terms of I = y/A like ro and r^ in the foregoing section. 
Solving the system of equations given by (|89p we obtain 
|o ~ —iy/l — E 2 , £5 ss iy/1 — E 2 and, therefore, £0 and 
£5 are complex conjugate. 

For the Taylor expansion of Ag we obtain the same 
formula (I5H|) with r replaced by £. However, in this case 
we get 



k 2 — k 2 — k 2 — fl 



for p = , 



(91) 



which means that the Lauricella function Fd reduces to 
an elementary function in the limit p = 0. For the linear 
correction with respect to a we get 



Ao, a (P) 



2nLE 



(92) 



The corrections with respect to n and I vanish, A$ >n (0) = 
= Agj(0), as well as the correction for A due to 
AeAO) = 0{l 2 ). 

D. Linear correction to frequency of 99 motion 



e e (0 = a 4 A(^ - 6)(e - 6) K - 1 1 U - 4 



In this section we will calculate the Taylor expansions 
of the integrals I vr and I v q defined in (J58|) , 




(89) 



where 0£(£) is the right hand side of ([5]) with the sub- 
stitution £ = cos 9, 



= (1 + a 2 A£ 2 - 4Aa<)(C - (n - a£) 2 )(l " £ 2 ) 
-(L-( a (l-£ 2 ) + 2<)£) 2 , (90) 



V-(p) = 2 



a § r (r)dr 



^«G0 = 2 



<$>e{0)d9 



The conversion of these two integrals to the standard 
form is more involved than in the case of A r and Ag 
because an additional function $ r or $g appears in the 



P{r) 



. The 



A,. 

0, which 



integrand here. Let us first consider $ r (r) = 

poles of $ r are given by the horizons A r 

we denote by hi, 1 < i < 4. Let /12 ~ and h 3 ~ 

2 correspond to the Schwarzschild case, whereas hi ps 

— Z _1 — 1 and /i 4 ~ / _1 — 1 tend to infinity for vanishing 

A. Then $ r can be rewritten as 



10 



fc P (r) 



aP{r) 

I 2 (r - ¥) (r - h 2 )(r - h 3 ) (r - ~^ 



(93) 



where hi = ?/ii, /14 = ^4. Now $ r can be decomposed in partial fractions and the resulting integrals transformed to 
the standard form by the substitution outlined in section IV Al This yields an integral of the form 



fcO^rl, (J. 
4 

4^ a/h 



(A r t + _B r )dt 



(1 - A^ r i) V*(i - *)(i - fepKi - fepKi - fep) 



Ar-F D f -,f3 2 ,2,m iipr \ + B r irF D I -,/3 2 ,l,m 



zyr 



(94) 
(95) 



where c^ and A^ vr are some constants with ATo vr = 0, ^2 = (5, |, |j l), and m^,, = (fc 2 r , fc| r , fc| r , Ni vr ). The Taylor 
expansion of this can be done analogously to section IV Bl and gives 



I V r(p) 



4Er 04 Il{n 2 ,k) 



(r 4 ~ 2Vr 03 (r 04 -r 02 )(l-£ 2 ) 



(96) 



as the other corrections with respect to n, I, and A vanish. Here roi, k, and n 2 correspond to the reference orbit, see 
section IPVl for a definition. 

Now let us consider I v g . With the substitution £ = cos 9 the integral can be transformed to 



I v e{p) = 2 



£3 



(L-(a(l-£ 2 ) + 2<)£)d£ 



s c 2 



(? - l)(e + 1) (C - i(2n* + /(n, /))) (£ - ±(2nZ - f(n, I))) ^/e^j 



(97) 



where f(n, I) — V4/ 2 n 2 — 1. Then the standard form can be obtained with the substitution described in section IV Al 
giving the integral (|95|) with r replaced by £. Here again c^ and N^ are constants with Nq v £ — 0. The Taylor 
expansion is given by 



I<p8 (p) ~ 2-7T signL — a , 



C 



(98) 



with vanishing corrections due to n, I, and A. 



E. Linear correction to frequency of t motion 

The two integrals to be considered for determining the correction to the t motion are 



, VW) 



e min VW) 



The procedure will be analogous to the foregoing subsection. The poles of T r and Tg are the same as of $ r and $0 
respectively and, therefore, the integrals can be expressed similar to (j9"5|) . We get 



Itrip) = 2 ^ Cj 



(=0 



A r -F D f ~,/?2,2,"W ) + B r nF D f -,/3 2 , l,m itr 



(99) 



/*r(0) + /tr, o (0)a + Jtr,n(6)n + W(0) A 



where Cjt r and A^ ir are constants with No tr = 0, f3 2 — (i, |, |, 1), mjt r = (fc 2 r , fc 2 r, ^Iri ^itr)> an d ^tr(O) = ^tr,o as for 
the reference orbit. The linear correction due to a is given by 



Wo) 



4£r 2 3 



(»'03('"04-ro2)(l -£ ;2 )) 2 (r 4-^03) 2 (' , 03 ~ r 02 ) 



r 4ViE 2 K(k) 2C(8 + 2C - 3C£ 2 ) 



(r 4-2)(r 03 -2) (r 03 - r 02 )(l - S 2 ) : 



25(fc) 



, (100) 
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where r\\ is a non-symmetric function of ro2, ?"03, and r 04i 

m = 2r 04 r 02 - r 03 rl 4 ro2 - 2r 2 4 r 03 + 2r§ 3 r§4 - 2ro 3 r 04 - ro 3 ro 4 r 02 + 2ro 3 r 02 
The correction due to n vanishes, /*,.,„ (0) = 0, and the correction due to A is given by 
E 

/tr,A(0) = 



V r 03(r04 - r 02 )(l - E 2 ) 



x < - ro 3 r 04 K(k) 
. r o 3 (ro4 -r 02 ) 



(o 2 + 12Cq 3 ) 



ai 



(1-25*) 




r ° 4 TT7 


,k) 


^02 + r 03 - 


3r 04 



£(*) 



24(l-£ 2 )(ro4-ro 3 ) 2 (ro3-ro2)(ro4-ro2) ' 12 

(a 4 + 12Ca 5 ) 8 - 75£ 2 + 149£ 4 - 116£ 6 + 32£ 8 



-i (r M -2)- 1 (ro8-2)- 1 



2AS 
238 - 754£ 2 + 952E 4 



e 2 {\-e 2 ) 2 

560£ 6 + 2 7 E 8 2C(2E 2 - 5) 



4C 
3 



(101) 



1-E 



2Y.I 



«6 + 4CQ7 



(l-£ 2 ) 3 



l-£ 2 



(l-£ 2 ) 2 



r Q3 r Q4: K(k) - r Q3 (r Q4: - r 02 )E(k) _ 2 II(ni, fc) 



'ro4 - r 3) 2 (r 3 - r 02 )(r 4 - r 02 ) 
256r 4 



r 04 -2 



n(n 2 ,fc) 



(102) 



where 5 = (r O4 -n)3) 2 (n33-n)2) 2 0""o4-n)2) 2 (n)4-2)(ro3- 
2) (^02 — 2) and o», i = 1, . . . ,7, are non-symmetric func- 
tions of ro2, ^03) and ro4 given in appendix [Bl 

The Taylor expansion for the 9 dependent part is much 
easier. The zero order term and all linear correction ex- 
cept I t e,a vanish, l t e(0) = lte,n(0) = Ite,A0) = 0. The 
remaining term is given by 



Ia,a(P) = 



2ttL 



c 



(103) 



I 

fundamental frequencies we will reduce the hyperelliptic 
integrals to elementary expressions and elliptic integrals, 
which can be handled easily with computer algebra sys- 
tems like Mathematica or Maple. In addition, we will 
perform the post-Newtonian expression for comparison 
with other results and use the parameter values of Mer- 
cury to get an idea of the order of magnitude of the cor- 
rections. 



Periastron shift 



VI. OBSERVABLES 



In section irni we defined the general expressions for the 
fundamental frequencies Q r , fig, and 0,^ in Plebariski- 
Demiahski space-time. The analytic expressions of these 
quantities are given in terms of complete hyperelliptic in- 
tegrals on a Riemann surface of genus 2 (see Eqs. (|76|). 
(|55|) , (|55j) , and (|99p ) , whose numerical evaluation is quite 
cumbersome. By a post-Schwarzschild expansion of the 

I 



In section ITVl we calculated the expression for fip for 
the Schwarzschild case, see (jT7|) . and noted that as long 
as a particle is considered to move in an orbital plane 
usually the difference angle Ap, see (|4"51) . is used instead 
of J7p. Let us now consider the linear corrections of $7p 
and Ap due to a, n, and A. As all linear corrections 
due to n so far vanished, the linear correction due to 
this parameter also vanishes for J7p and Ap. The linear 
correction due to a for the frequency ilp is given by 



n ,=*, W0) + sign! VCA r . a (6) - EA r . (l - M) + j ( ^ )(27r _ y/CAr,o) 
fip,a(0) = + signL ; 



itrfi 



T 2 

1 tr,0 



r 4 „, , , , /^ r 1 Ci(^Q2 ~ r 03 )(g 2 - l) 2 A-(fc) - 2r 03 (4 - 3C(E 2 - l))E(k) ( \L\ 

-n(n 2 , k) + VC\L\ -y- -2-— ^ — -T3 K(k) 1 - -= 

ro4 - 2 (r 02 - r 3) z (r 2 - ro4j(ro3 - r 0i y(E* - If \ y/C 



\L\ 



E 2 Z 2 



K(k) 



r 03 r 4mE 2 K(k) 



2C(8 + 2C-3CE 2 )E(k)\ 



(r 4-2)(ro3-2) +r ° 3 (r 03 - r 02 )(l - E 2 ) 2 ) 



x ((1 - E 2 )(r 04 - ro 2 )(r 04 - r 03 ) 2 (r 3 - r 02 )) 



2TrVro3{ro2-r 04 ){E 2 -l) - iVCK(k) 



(104) 



where Z = -r 03 r 04: K(k) + r 03 (r 04: - r 02 )E(k) + 2r 4 3 1 _ 2 ^ n(m,fc) + ^^n(n 2 , k) and Ci, »7i are non-symmetric 
functions of ro 2 , ^03, and ro4 defined in (|101[) and (1831) . Here r"oj, fe, and rij correspond to the reference orbit and 



12 



are defined in section HVl Note that Op ja does not depend on the sign of L, i.e. on whether the particle travels on a 
prograde or retrograde orbit. 

The correction of Ap due to the angular momentum of the gravitating source reads 



A P ,„(0) = 



4EVC 



y/ '(1 - E 2 )r 03 (r 0i - r 02 ) 



\L\Ci 



{ro4 - r 3) 2 (r 03 - r 02 )(r 04 - r 02 )(E 2 - 1) y/C 



2|L|r 03 (4-3C-3C*£ 2 ) 



(1 - E 2 ) 3 (r 0i - r 03 ) 2 (r 03 - r 02 ) 2 (r 4 - r 02 ) 



E(k) 



7"04 



ro4 



-n(n 2l fc) 



c 



K(k) 



(105) 



r 



which also does not depend on the sign of L. These for- 
mulas are exact in M and, therefore, quite complicated. 
If we consider the post-Newtonian approximation of these 
terms by using ([2? ]) -([30 ]) and assuming that r 03 and r 4, 
the turning points of the reference orbit in Schwarzschild 
space-time, become large we obtain 



O P , a (0)a = 2 



1 — 6 sin i 



d 3 (l 



^JM 2 



A P ,a(0)a = 47T 



1 — 6 sin #01 
di(l-e 2 )l 



JM* 



0(JM 3 ) , (106) 
+ 0(JMi), (107) 



where d is again the semimajor axis, e the eccentricity, 
and J = aM is the non-normalized specific angular mo- 
mentum of the gravitating source. These terms are sim- 
ilar to the expression derived by Lense and Thirring [4] 
for the precession of the longitude of pericenter, to which 
it differs only in the second term containing the inclina- 
tion #oi • However, the fundamental frequency Op is not 
identical to the precession of longitude of the pericenter: 
As O r does not contribute to the first order correction 
of Op the value of Op a is calculated by the correction 
due to O v — -j£ only, which is the averaged -£. The 
longitude of pericenter on the other hand contains the 
argument of the pericenter which is defined not on the 
reference plane but the orbital plane. 

I 



Note that the expressions for Op and Ap depend on 
the inclination of the orbital plane. For an inclination of 
#oi < arcsin i w 0.0537T the periastron shift is perturbed 
in the direction of the rotation of the gravitating source 
whereas for <?oi > arcsin g the perturbation is in the op- 
posite direction. (Here an inclination of 0qi = § corre- 
sponds to the equatorial plane and 6>oi < ■§ by definition, 
see sec. HVl ) In particular, if the unperturbed test parti- 
cle moves on the equatorial plane, the periastron shift is 
perturbed against the direction of rotation. This seems 
to be counter-intuitive as a particle radially approaching 
the gravitational source is dragged along the direction of 
the rotation. However, this does not mean that the shape 
of the orbit is effected in the same way: the rotation of 
the gravitating source acts as an repulsive force which 
also causes the peri- and apoapsis to increase, see Q. 

Let us now consider the linear correction of Op and Ap 
due to A. For Op we get a very complicated expression 
given in appendix [Bj In terms of the expression derived 
in the previous section we get 



Op,A = signL v C 



A 



r.A 



'tr,0 



27T \ Itr.A 



(108) 



The correction of Ap due to A is much simpler, 



A P ,A = 



signL \[C J 2C(r 02 ro 3 + r 02 ro4 - 2r 03 ro4) - 71 

v /(l- J B 2 )ro 3 (ro4-ro2)(l - E 2 ) \ (r 4 - r 03 ) 2 (r 4 - r 02 )(ro3 - r 02 ) 



ro3r 4K(k) 



48 + 16C - 76CE 2 + 57CE 4 + C 2 - C 2 E 4 
16 + 8C - MCE 2 + 27CE 4 + C 2 - C 2 E 2 



ro3(»"o4 - r a2 )E(k) 



6r 04 
1-E 2 



n(m,fc) 



(109) 



These complicated expressions can be simplified by con- 
sidering the post-Newtonian approximation which reads 



A P , A A = 7rsignLtfVl-e 2 77 



2 - e 2 - 
9.TT*\ Z nT,d 2 A +rO(AM), (111) 

\/l — e 2 



fi P , A A= isignLdVl-e 2 AAf5 +0(AM%), (110) 



where A : 



3A 



jgz is the usual non-normalized cosmological 
constant. The first term of Eq. (|111[) coincides with the 
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result in [33| . Interestingly, although we assumed both A 
and M to be small we assumed nothing about their ratio. 
This raises the question whether the first term of (jllip 
is indeed the largest of the expansion. But if we divide 
the second term by the first this is proportional to -t-, 
which is assumed to be small. This proves that the terms 
of the above expansions in fact decrease in magnitude for 
increasing powers of M. 

The exact formula for the perihelion shift assuming 
a non-vanishing cosmological constant was derived in 
[Ull. We will use the formalism of P3 (see Eq. (70)) to 
compute the perihelion shift of Mercury in Schwarzschild- 
de Sitter space-times and compare it with our approxi- 
mate formula. For this, we use M = 1476.625 28 m and, 
by using (|2"T)) - ([3"0")) . we determine averaged values for the 
energy and angular momentum from the apo- and peri- 
apsis data given by [3CJ] and written down in IIVE1 With 
the exact formulas given in [12j the difference between 
the perihelion shift for A = and with a cosmological 
constant of A = 3A/A/ 2 = 3 x 10~ 52 m~ 2 is given by 



Ap,(A=io- 52 )/yM — Ap,(A=o)/y 



M 



1.038 833 075 425 928 x 10" 



arcsec 



cy 



(112) 



where y M is the Mercury year as calculated in (|56|) . If 
we insert the same values in our formula (|109[) we get 



Ap/y M - Ap,o/y 



M 



= 1.038 832 108 177 831 x 10" 



cy 



C(A 2 ) (113) 



which agrees very well with (|112[) . 



B. Lense-Thirring effect 

The Lense-Thirring effect, which can be identified with 
a precession of the orbital plane in the weak field limit, is 
up to first order not influenced by any other parameter 
than the rotation a. This is because all linear correc- 
tions to Tg and T v due to n and A vanish. The linear 
correction due to a is given by 



n 



LT> = 



2r 04 



(ro4 



-U(n 2 ,k)-2K(k) 



(114) 



where Z = -r 03 r 04: K(k) + r 03 (r 4 - r 02 )E(k) + 

2r 04 tzw n ("-i; fc ) + ^7^2 n ( n 2, k) andr 0l , k, and rn cor- 
respond to the reference orbit, see section HVl 

Let us consider now the weak field limit by assuming 
that r 03 and r 4 become large. Using (l27l) - (j3"0"|) we get 



^LT, Q a — 



d 3 (l 



JM 2 + 0(JM 3 ) . 



(115) 



where J = aM is the non- normalized angular momentum 
of the gravitating source per unit mass, d is the semima- 
jor axis, and e the eccentricity. This formula is identical 
to the precession rate of the longitude of the ascending 
node as given by Lense and Thirring, cf. [4] . 



C. Conicity 

The effect of n can be observed by a deviation from 
the symmetry of the geodesic motion with respect to 
equatorial plane. In sec. IIIII we defined A con i C i ty = 
7T — (9min + ma x) as a measure for this deviation, which 
is given by 



^conicity 



4E 



(116) 



up to first order. Linear corrections due to a and A van- 
ish. This deviation from the symmetry to the equatorial 
plane also implies that for n/Oa test particle can not 
move on an orbital plane and, in particular, not in the 
equatorial plane (if £^0). Instead it moves on a cone 
with opening angle n — -^j=n + 0(n 2 ). As 6> i and 8 02 

are perturbed by the same value a = — signL -^n (up 
to first order) the symmetry axis of this cone coincides 
with one of the two normals of the orbital plane in the 
Schwarzschild case, i.e. 8 n = 9 n2 — § G (0, -|) if a < or 
0n = 0oi + f e(f,7r)ifa>0. 

In the post-Newtonian limit the conicity becomes 



A 



conicity 



I signL 



4 signL h 
^d{l~e 2 ) Mi" dl(l-e 2 )i 



nM5 , 



(117) 



where h = nM is the non-normalized NUT parameter. 
Here d and e are the semimajor axis and the eccentricity, 
as before. Note that as we assumed n = -tj to be small 
in the first place, the first term in the above equation will 
indeed tend to zero for M — > (but is still the largest of 
the expansion). 

The conicity can be used to determine an upper bound 
for the value of n. As the error tolerance of the inclination 
of an orbiting object places a bound for A con i c it y , this can 
be used in one of the above equations to get an estimate 
of n. From [30] we obtain for the difference of maximal 
and minimal inclination of Mercury over a time span of 
10 years | A conicit y| < 4.200 1 arcsec or | A conicit y| < 2.04 x 
10" 5 rad. Inserted in (flTBl this yields |n| < 0.032. 



VII. CONCLUSIONS AND OUTLOOK 

In this paper we derived analytic expressions for the 
observables periastron shift and Lense-Thirring effect as- 
suming neutral test particles in the general axially sym- 
metric Plebahski-Demiahski electrovac space-time with 
vanishing acceleration of the gravitating source. We de- 
termined the direct dependence of these observables on 
the parameters of the space-time by expanding them in a 
Taylor series up to first order. In addition, we defined the 
conicity of an orbit and analyzed its dependence on the 
space-time parameters by also expanding it in a Taylor 
series up to first order. Within the linear approximation 
we used the osculating orbital elements of Mercury to 
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derive an upper bound for the dimensionless NUT pa- 
rameter n. 

From the six parameters of a Plebahski-Demiahski 
space-time with vanishing acceleration, only the Mass 
M, the rotation a, the NUT charge n, and the cosmo- 
logical constant A have any linear effects on neutral test 
particles. For Schwarzschild space-time (0 = a = n = 
A = Q c = Q m ), the only effect is the periastron shift, 
whereas the Lense-Thirring effect vanishes and the mean 
value of the polar coordinate coincides with the equator. 
In the following we summarize the effects due to the other 
parameters compared to Schwarzschild space-time: 

• Taub-NUT space-time (0 = a = A = Q c = Q m ): 
The mean value of the polar motion deviates from 
the equator and the motion takes place on a cone 
rather than a plane. 

• Kerr space-time (0 = n = A = Q c = Q m ): The 
periastron shift is changed and the Lense-Thirring 
effect is nonzero. The latter can be interpreted as 
a precession of the orbital plane in the weak held 
approximation. Both are independent from the di- 
rection of rotation of the particle. 

• Kerr- Taub-NUT space-time (0 = A = Q c = Q m ): 
In addition to Kerr the mean value of the polar mo- 
tion deviates from the equator. The combination 
with a non-vanishing Lense-Thirring effect causes 
a precession of the orbital cone. 

• Kerr-Taub-NUT-de Sitter space-time (0 = Q c = 
Q m ): In addition to Kerr- Taub-NUT the periastron 
shift is changed by a non- vanishing A. However, 
due to the smallness of A the effect is tiny. 

The effects of non- vanishing space-time parameters may 
also be analyzed compared to any other space-time cov- 
ered by the general Plebahski-Demiahski metric using the 
same methods as presented in this paper. In particular 
the Kerr space-time may be used as reference. In this 
case Tg y0 and T Vy0 are also given in terms of elliptic in- 
tegrals instead of elementary expressions. 

The mathematical framework presented in this paper 
is based on the theory of hyperelliptic and elliptic func- 
tions. We defined a standard form for hyperelliptic inte- 
grals based on the Riemann form of elliptic integrals and 
expressed it in terms of the hypergeometric Fd function 
introduced by Lauricella (see also [26|]). The methods 



presented in this paper may also be used to analyze the 
linear (and higher order) effects of the six parameters 
of Plebahski-Demiahski space-time on charged particles 
with a complete analogous procedure. Another inter- 
esting question is whether additional effects arise for a 
higher order approximation, for example the quadratic 
effect of n on r-dependent expressions or the coupling 
of a to the other parameters. Similarly to the analysis 
presented here observables for unbound orbits like light 
deflection or the deflection of massive particles may be 
considered. This will be postponed to a later publication. 
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Appendix A: Lauricella's Fr> function 

The four functions Fa, Fb, Fc, and Fd of Lauricella 
are hypergeometric functions of multiple variables gen- 
eralizing the hypergeometric functions of Gauss and Ap- 
pell. They were introduced in 1893 j35| and given as a 
hypergeometric series 

F D (aJ, 7 ,x) = £ ^L% x*, (Al) 

where t*is a multi-index, \x Ln \ < 1 for all n, and (•) the 
Pochhammer symbol. Here |t| = J^ n i n , V. = J\ n t n !, and 
{fi)v = Y\ n (Pn)i n - The function Fd can be extended to 
other values of x by analytic continuation. 

In this paper the Fd function is used because it can 
be represented in an integral form 



F D {a,/3,y,x) 



r(7) 



r(a)r(7 — a) 
x [ i"- 1 (l-i)''- Q - i n(l-.x„i)- /3 "di (A2) 

for Re(7) > Re(a) > 0, which is exactly the form of all 
hyperelliptic integrals appearing in this paper. 
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Appendix B: Details of fip, 



The calculation of the linear correction of the periastron shift due to the cosmological constant yields a very 
cumbersome expression given by 



n 



P,A 



2£signLVC 
[E 2 - \)Z 



71 - 2C(r 2r 03 + r 02 r 4 - 2ro 3 r 04 ) 



(V 2 - r 04 )(r m - r 03 )(r Q3 - r 4) 2 



ro 3 r iK(k) 



-16C + 76CE 2 - 57C£ 4 - C 2 + C 2 E 4 - 48 
16 + 8C - 36CE 2 + 27CE 4 + C 2 - C 2 E 2 



ro3(ro4 - r 02 )E(k) 



6r 04 
E 2 -\ 



n(m,fc) 



signL 



2VCK(k) - TTy/r 03 (ro2-r M )(E 2 -l) 



2EZ 2 

j r 02 + r 03 - 3r 04 
J (£ 2 - l) 2 

2 



r 03 r 4-ftr(fc) - r 3(r 4 - r m )E(k) - 2r Q4 _ 2 II(ni,fc) 



ai 



L(£ 2 -l) 3 12(r 4-2)(r 03 -2) 

a2 + 12Ca: 3 



+ 



24(r 04 - 2)(r 03 - 2)(r 02 - r 04 )(r 2 - r 03 )(r 03 - r 04 ) 2 (E 2 - 1) 
4C a 4 + 12Ca 5 



r-03^04^(fc) 



3(£ 2 _ 1)2 24(r 04 - 2)(r 03 - 2)(r 02 - 2)(r 02 - r 04 ) 2 (r 2 - r 03 ) 2 (r 03 - r 04 ) 2 



-75£ 2 



149£ 4 - 116£ 6 + 32£ 8 



E 2 (E 2 _ !)4 
2C(2£ 2 _ 5) at , _|_ 4Cq, 7 



^03(^04 -ro 2 )E(k) + 



256ro4 
ro4 - 2 



II(n 2 ,fc) 



(£ 2 - l) 2 8(r 02 - r oi )(ro2 - r 03 )(r a3 - r 0i ) 2 (E 2 - 1) 

238 + 128£ 8 + 952E 4 - 754E 2 - 560E 6 ' 



(E 2 1) 



r 04 n(ni,fc) 



(Bl) 



where Z = — ro 3 ro4-ff (fc)+ r 03( r 04~ ro2)E(k) + 2ro 4 3 1 J^ II(m, k)+ r Sr °^ 2 H(n 2 , k) and 71, ai, . . . , aj are non-symmetric 
expressions of r 02 , r 03 , and r 4. These constants read 



2 „2 



„2 „2 



71 = -2r 03 r 04 + 3r 03 r 4 + 3r 04 r 03 - 3r 03 r 02 - 3r 04 r 02 + r 03 r 02 - 2r 03 r 04 r 02 + r 04 r 03 r 2 + r 04 r 02 + r 0Z r M r 02 (B2) 



ai =1536 - 288r 04 - 288r 03 + 96r 02 + 20r 2 2 - 8rQ 4 r 02 + 4r 2 3 r 04 r 02 + 4r 2 )4 r 03 r 2 + 5r 03 r 04 r 2 2 

+ 8r 02 r 03 r 4 + 14r 2 3 r„ 4 + 15rg 4 r 03 - I0r 03 rl 2 - I0r 04 rl 2 - 30rg 3 - 8r 2 3 r 02 + l^rl 4 r 03 + 14r 2 3 r 4 



+ 56r 03 r 04 - 32r 02 r 04 - 32r 02 r 03 - 84r 2 3 + 15r^ 3 r 04 - 30r, 



04 



84r, 



04 ; 



(B3) 



a 2 =864r 03 r 2 , 4 r 2 2 - rf )3 r 4rl 2 - 336r 03 rj] 4 ro 2 + 384ro 3 r 04 ro 2 - 20rl 3 r 04 rf )2 + 384rl 4 r 03 r 2 + 36ro 3 ro4ro 2 
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■' " :1 '"■■■-'■■-■■ -->.'...;.. 290r^r 2 ,r 2 o +432r3 3 r 04 r 02 + 23^ 3 r^ 2 r 4 
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76r 
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816ro 4 r 02 

(B4) 



16 

a 3 =32r 2 l4 r 02 + 2r 03 r 2 4 r 2 2 + 14r 03 ro 4 r 2 - 32rl 3 r 04 r 2 - I6r 03 r 04 r 2 2 + 4r 2 3 r 02 ro 4 + 3r 03 rl 2 rl 4 - 4r^ 3 r 2 4 r 02 
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- 168r 2 3 r6 2 - 168r 6 2 r 2 4 + 32r 6 3 r 4 4 + 204r 7 3 r 2 2 + 640r^ 2 r 4 4 - 220^ 2 r 3 4 - 340r^ 4 r^ 2 + 584r 6 4 r 4 2 - 414r 7 4 r 3 2 

- 174r 7 2 r 3 4 + 152r 4 4 r 6 2 + 128r 5 3 r 3 2 + 348r 2 3 r 7 2 + 672r 4 2 r 3 4 + 1664r 3 3 r 5 4 - 1896^ - 268^^ - 102r 7 3 r 3 2 
+ 80r 4 3 r^ 2 + 736r^ 3 r 4 4 - 832r 5 3 r 3 4 - 264r 6 3 r 2 4 - 592r 4 3 r 4 4 + 4608r 2 2 r 2 l3 r 2 l4 + 1536r 3 2 r 3 3 - 96r 3 3 r 4 2 

+ 864r 5 3 r 2 2 + 1632r 5 4 r 2 2 - 768r 4 2 r 2 4 + 1536r 3 3 r 3 4 - 2400r 4 4 r 3 2 - 864r 4 3 r 3 2 , (B7) 
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a 5 = - 18r^ 2 r^ 4 + I28r 03 rl 4 rl 2 - 16r 4 3 r 4ro 2 - 128ro 3 ro 4 r- 2 - ^Sr^ 2 r Q4 r Q3 + 24^3^4^ + 24^3^4^ + 32ro 2 rg 4 ro 3 
+ 24r 4 3 r 04 r ! 2 - 40r^ 2 r 4^ 3 - 8r^ 2 r^ 4 r 03 + 56r 03 rj] 2 rj] 4 + 16r^ 3 r^ 2 r^ 4 + 120r^ 3 ro 2 rg 4 - 40rg 2 r^ 4 r 03 

- 168r 03 ro 2 ro4 - 56r 4 4 r 2ro 3 - 88r 03 ro 4 ro2 + 24r 03 r 4 4 ro 2 - 24r^ 3 ro 4 ro2 - 24r 4 3 ro2ro 4 + !20rQ 3 r 2 4 r 02 

- 32ro 3 ro 4 r 02 - 16^3^4 r\ 2 - 128ro 3 r 04 r 02 + ^l 3 rl 2 r M + Qrl 3 rl A r 02 + 176r 03 r Q2 rf )4 + 32rg 2 r Q4 r 03 + 48r 4 3 7- 4ro2 



128r 2 3 r 2 4 r 02 + 128r 2 3 r 2 2 r 04 - 18rg 3 rg 2 - W4rl 3 r 04 r 2 2 + 20r 4 2 r 04 r : 3 



247'o2 r 04 



32r 03 r 04 r i 



+ 36^3^2 - 24r 4 3 r ! 4 + 22r 02 r^ 4 - 12r 02 r 3 03 ri 4 + 20r 4 02 r 3 04 r 03 + 20r 4 02 r 2 4 r 2 3 + 32r { 



3 „2 3 

03' 02' 04 



02 

64r, ; 



^ r Q3 r Qi r Q2 

3 r 2 
02 '04 



48r Vb 3 2 ro ! 4 + 96r^ 2 rg 4 + 16^ 2 ^ 3 r^ 4 + 22ro 3 rg 4 r§ 2 + 4r 4 3 r : 4 ro2 



88r^ 2 r 4 4 



88r&rS4 



+ 64<or^ + 12r, 



3 '04 



03 r 04 + 4r 03 r 04 + 36r 02 



24^,^ - 8r 4 2 rg 3 - 64rj^ 2 + 36^ 3 r 4 4 + 64r 3 l3 r ! 4 + 3r{j 3 ro4rg a 



44 ^o. 

11^05^04^02 + GroVoW + 8r 4 3 r(S 2 - 6r^r^ 2 ^ 4 - 3r^ 2 rg 3 r 4 4 + 9r^ 2 rg 4 r 3 - 12r£ 2 rj^ 3 - 4r 4 3 r 4r 4 2 



+ 6r 03 r 04 r 02 + 9rg 3 r 04 ^ 2 - 7r 02 r 04 rg 3 + 3r 03 r 04 r 02 - 2r 02 r 4 3 r 04 + llr 02 rg 3 r 04 



0^4 2 3 

O '0 3 '02'04 



r- 4 r- 3 r 2 
'02'04'0 3 



2^n 3 r 02 r n4 + 9 r o 3 r 04 r 02 ~ 12?"o 3 r n4 r n2 + 1 ' r 02 r o 3 r o 



'03' 04' 02 



02'0 3 '04 



-6r^ 4 + 4r 4 3 r 4 4 -48rj 

" ^ 4 03 r 3 02 + 64r 3 04 r 2 2 + 64r 2 2 r 3 03 , 



dZrr.^r, 



02 '03 



'02'0 3 '04 
12 r 3 r 02 



22^02 r 3 r 04 



22r 04 r 02 



3 

03' 02 



QrLr, 



44r 4 2 r 3 14 



22?*o 3 r 04 



12r, 



3 '02 



44»Vo2 - 8r 4 2 r 2 4 



-48r; 



3 r 04 



4rir, 3 



04' 02 

(B8) 



a 6 =32r 03 r ! 4 ro2 - 58r^ 2 r 04 ro 3 + 128ro 3 r 4r 2 - 28ro 3 r 4ro 2 



128ro 4 r 03 r 02 + 5rl 3 r 04 rl 2 + 26r 3 12 ro4 r o 3 



35r 4 3 ro 4 ro2 - 1 8 r^r 04?^ 



18?-02 r 04 r 3 



67r 03 ro 2 rg 4 - 64rl 2 r 04 r 03 



r 3 r 02 r 04 



256r 03 r 4r 02 

+ l 28r o 3 r 02?"o4 + 68r 02 r 04 r 03 - 260r cl3 r ; 2 r i 4 + 5r 4 4 ro 2 ro 3 + 49r 03 r^ 4 r 02 + 5> 03 r 4 4 rg 2 + 9rQ 3 r 4ro 2 
- 39r 4 3 ro 2 ro4 - 160rl 3 rl 4 r 02 + 30r^ 3 r^ 4 r 2 + §rl 3 rl 4 rl 2 + 32r 5 3 ro4? , o2 + 140r 03 ro 2 r 4 4 - 88r 4 2 r 04 ro 3 



29r2,r? 



IQrLr 2 



52r 03 r Q4 r Q2 - 4r 03 r 4r 02 -r *j» 03 , 02 -r iu. 03 / 04 
236r ' 4 r 2 + lMr 2 2 r 4 4 - 304^ 3 r^ 4 + 32r 4 4 r 03 



69r 04 r 02 - 17r 03 r 2 - 256r 03 r 04 



32?'o2 r 'o4 



- 36rLr, 



3 

02' 04 



327"o 3 r 04 



H r 3 r 04 



15?"o 3 r 04 



b27Y,Tr n4 



69rg 4 r 3 + 29r, 



.5 r 2 

02' 04 



36r, 



03 r 02 



6 r 02 r 04 



44r 4 2 r ! 3 



32r 3 3^ 4 + 128r 3 4 r 03 
I28r 2 4 r 2 2 + 2r% 3 r 2 2 + 128r 2 
48r 4 4 r 3 12 + 236r 03 rg 4 



128r3 3 r 02 
- 12r^ro2 + 12rg 3 r 4 - 
2 + 128r3 3 r 04 - 34r, 

12 r 3 r 02 7 



32r, 



2 „3 

3 '02 



32r 4 3 r 02 - 58r, 



3 '04 



128r 02 r 03 r 04 4r 02 / 03 



k>3 ™o 



5 „2 

04' 02 



32r 04 r 02 + 44rQ 2 r 04 



6r, 

f 168r 3 3 r 04 



3 '02 
2 1 icoJ „3 



- 128r3 4 r 02 
17r^ 3 r 4 
■ 32r 4 3 r 04 



(B9) 



r 



These complicated expressions are exact and may be sim- 
plified by an additional approximation. For example, for 
large peri- and apoapsis only the higher orders of ro 3 



and r$ 4 may be taken into account. In particular, in the 
post-Newtonian case these long expressions are reduced 
to (fTTUD . 
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